Spin-gap phase of a quantum spin system on a honeycomb lattice 
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We study a quantum spin system on a honeycomb lattice, when it includes frustration and dis- 
tortion in antiferromagnetic (AF) exchange interactions. We transform the spin system onto a 
nonlinear a model (NLSM) in a new way preserving the original spin degrees of freedom. Assisted 
by a renormalization-group argument, the NLSM provides a ground-state phase diagram consisting 
of an ordered AF phase and a disordered spin-gap phase. The spin-gap phase extends from a strong 
frustration regime to a strong distortion regime, showing that the disordered ground states are es- 
sentially the same in both the regimes. In the spin-half case, the spin-gap phase for the spin system 
on a honeycomb lattice is larger than that for the J1-J2 model on a square lattice. 
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Recently, materials regarded as spin systems on honey- 
comb lattices have been found. Kataev et al. 'I'l synthe- 
sized a spin-half {S — ^) material on a honeycomb lat- 
tice, InCu2/3Vi/303. They estimated the first-neighbor 
exchange parameter as 280 K, and found an antiferro- 
magnetic (AF) ordering at 38K with interlayer interac- 
tions. Their results are confirmed in refined samples by 
Kikuchi et al. 0. Miura et al. reported other mate- 
rials, NaaTaSbOg (T = Cu, Ni, and Co). The material 
with T = Cu is an 5 = i system and shows a spin-gap 
in magnetic susceptibility. The materials with T = Ni 
and Co are 5 = 1 and S — ^ systems, respectively, and 
both have AF order. So far, samples are reported to 
have fairly strong dimer-like distortion (T = Cu) 4] or 
randomness (T = Ni and Co) in exchange interactions on 
honeycomb lattices. In view of recent experimental activ- 
ity, materials regarded as various types of spin systems 
on honeycomb lattices are expected to be synthesized. 
Hence it is worthwhile to show a theoretical scheme for 
spin systems on honeycomb lattices. 

In two dimensions, there are two kinds of simple bi- 
partite lattices, a square lattice and a honeycomb lat- 
tice, if only first-neighbor AF interactions are considered. 
A spin system on a bipartite lattice has a ground state 
with an AF order in a large spin-magnitude limit. A 
question is whether or not a disordered spin-gap phase 
appears with destroyed AF order, when the spin mag- 
nitude becomes small so that quantum fluctuations are 
large. Quantum fluctuations are further enhanced by 
frustration among AF interactions and by distortion as 
inhomogeneity in interactions. Hence a disordered state 
will be interestingly studied under variable strengths of 
frustration and distortion. 

In the square-lattice case, the J1-J2 model has been 
eagerly studied. Theoretically, whether or not a disor- 
dered state is formed around J2/J1 = 0.5 is a problem 
controversial for long time. The difficulty originates from 
the fact that the region of the disordered phase is very 
narrow in the J1-J2 parameter space even if it exists. 
Further a material realizing the J1-J2 model in a disor- 
dered state has not been found. The difhcult synthesis of 
such a material also seems to come from the narrowness 
of the disordered phase. 



In contrast, a quantum spin system on a honeycomb 
lattice is more hopeful to find a disordered ground state. 
In fact, the coordination number 3 of a honeycomb lattice 
is smaller than 4 of a square lattice, and so the AF order 
for a honeycomb lattice is more fragile. The case without 
distortion is treated by a numerical diagonalization and a 
linear spin- wave theory ITll| , Schwinger-boson mean- field 
theory 10], and a nonlinear a model (NLSM) method |^. 
A systematic study of a spin system on a honeycomb 
lattice both with frustration and distortion is expected. 

Among various methods to analyze a two-dimensional 
spin system, mapping it onto an NLSM and applying a 
renormalization group (RG) analysis is effective to study 
a quantum phase transition. In the case of a square lat- 
tice, a spin system only with first-neighbor exchange in- 
teractions is shown to be mapped onto an NLSM without 
topological term [||. For a J1-J2 model, NLSM methods 
are developed and a spin-gap phase are studied Af- 
ter then, a refined mapping onto an NLSM is proposed 
for a J1-J2 model on a square lattice with plaquette dis- 
tortion [g. In the derivation, the original spin degrees 
of freedom are properly treated without doubling by a 
variable transformation, in which a slowly varying AF 
variable appears for a block of four spins. The derived 
NLSM has no topological term even if a second-neighbor 
interaction and a plaquette distortion exist. The RG 
analysis is successfully applied to the obtained NLSM 
and provides a phase diagram in the space of parameters 
describing frustration and plaquette distortion 

In this paper, we propose a novel NLSM method for a 
spin system on a honeycomb lattice in the case with frus- 
tration by second-neighbor AF interactions and dimer- 
like distortion in first-neighbor AF interactions. The 
original spin Hamiltonian is mapped onto a field theory 
in a manner that six close spins are transformed into one 
slow modulation of the AF spin configuration and five 
rapid motions. The degrees of freedom of the original 
spin model are properly considered through the deriva- 
tion, and the choice of the cutoff becomes unique . The 
derived field theory is an NLSM except for an additional 
term. After dealing with the additional term, we apply 
the RG analysis Qf to the NLSM to obtain a ground-state 
phase diagram; We obtain a ground-state phase diagram 
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FIG. 1: Spin system on a honeycomb lattice, where a spin 
with magnitude S is on a lattice point. A first-neighbor 
exchange interaction with coupling constant Ji (J[) is rep- 
resented by thin (bold) solid line. A second- neighbor ex- 
change interaction with coupling constant J2 is represented 
by a dotted hne; dotted lines should be between any second- 
neighbor spins, although not drawn to avoid complication. 
The pth block is shown by a dash-dotted line; the adjacent 
blocks are labeled as p„ {j = l,2,3;u — l,r). Variables n,, 
{j = l,2,3;u = l,r) represent spins in the path-integral for- 
mula with spin coherent states. 

in the parameter space, and show a disordered region in 
which there is a spin-gap. 

The Hamiltonian of the present spin system on a hon- 
eycomb lattice is 

i/ = ^J.,S,.S„ (1) 

where Si is the quantum spin at site i and the spin mag- 
nitude S is uniform and arbitrary. Exchange constant 
Jij takes Ji, J[ or J2 depending on the bond between 
sites i and j as shown in Fig. ^ The distance between 
first-neighbor spins will be written as a. We consider 
the quantum Hamiltonian in the classical Neel ordered 
region. 

Quantum fluctuation is increased by distortion and 
frustration in exchange constants. To measure the 
strengths, we introduce distortion parameter i5 and frus- 
tration parameter a as follows: 

Ji = Ji(l - 5), J[ = Ji(l + 2,5), a = J2/J1. (2) 

Then Ji is an average of all first-neighbor exchange con- 
stants: Ji = i(2Ji + J[). 

We construct a path-integral formula by means of spin 
coherent states. For a spin coherent state labeled by a 
unit vector Hjir) at imaginary time r, the expectation 
value of Sj is represented as (Sj) = CjSnj{T), where Q 



takes -|- or — depending on the sublattice which the jth 
site belongs to. Then the general formula of the partition 
function at temperature is written as 

Z=jD[n,{T)] Y[5{{n,{T)f-l)e-^ (3) 
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with action 

A^iSY,(:M^A+ / dTH{T). (4) 

The first term in A is the Berry phase term with the solid 
angle w[nj] which the unit vector nj(r) forms in period 
/?. Hamiltonian H{t) is written as 

H{r) = \s^Y.J^>^^^)-^^^^)f^ (5) 

using constraint (nj(r))^ — 1. Hereafter we do not ex- 
plicitly denote the r dependence of nj(T). 

By finding an appropriate variable transformation, we 
will separate out a slow mode describing AF fluctuation 
for the honeycomb lattice. We adopt the region enclosed 
by a dash-dotted line in Fig. as a unit of transforma- 
tion, and call it a block; the first-neighbor blocks to the 
pth block are labeled by p{ (j = 1, 2, 3 and u = rj). 
Accordingly we relabel the six variables, rij's, in the pth 
block as ril{p) and n;^ (p) {j — 1, 2, 3), as shown in Fig. 
^ Then we perform the following variable transforma- 
tion for each block: 
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Here L(p), X±(p), and Y±{p) describe fluctuation modes 
around AF variable ni(p). Six original constraints, 
{nl{p)y = 1 (j = 1, 2, 3 and t = r, I), are changed 
to six new constraints, (m(p))^=l and m(p)-L(p) = 0, 
m(p)-X±(p) = 0, and m{p)-Y±{p) = [H. 

Taking a continuum limit, the first term of action Q 
is written as 

p i=i 
2 Q f 

^ -i — dTd^r(3L + X+)-(mxdrm). (7) 
3a J 

In each term of Eq. ||SJ|, calculation like the following is 
carried out; 

nlip) - nfipl) = m(p) - ni{pl) + 6a(L + X+) 
^ 3a[5:rm-|-2(L + X+)] (8) 
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up to the lowest order of derivatives and fluctuation vari- 
ables. Thus the lattice action becomes to the field- 
theoretic action: 

A = [[(dTd^v[-^{J[-U2){d,mf 



+ B 

2X, 



+ 



(B^ 



-4J2)(a,m)^ 

+ 18 JiL^ 6L • (2 JiX 

+ 2( Ji -t- J{ - 3J2)X^ 

+ 2(Ji - 3J2)(3X^ - V3X_ • dyxn 

+ 6(J(-3J2)Y2} (9) 

with B+ = i[3J( -t- i(aS')-im X 9^m], B_ ^ 
\[{Ji - J() a^m - i(aS')-^ mxa^m], Ji = i(2Ji-}-J{), 
and J{ = \{Ji + 2J(). This action includes all the low- 
energy excitations surviving the continuum approxima- 
tion, since the original degrees of freedom are not spoiled 
in the variable transformation ((HJ. 

Integrating out the partition function with the action 
with respect to massive fields L, X±, and Y±, we 
obtain the following action: 



A' = 



+ 



1 



3a 



-{dyUl) 



(10) 



with c+ = 1 + (5 — 2(5^ — 6a and c_ = 1 — 5 — 6a. Except 
for the first term this is an NLSM action. The first term 
is rewritten as 

.SS 



iAr 



3a 



dyQiy) 



with 



Q{y) — jj drdx m • {drTxixdxVcv)- 



(11) 



(12) 



The quantity Q{y) is not a topological charge in the 
three-dimensional Euclidean space. For a fixed value of 
y, however, Q{y) is the topological charge of the field m 
in the Euclidean t-x plane, and is an integer. Further, 
for the continuity, Q{y) must be the same integer for all 
values of y. This means that there is at least a vortex 
line for Q 7^ 0. Hence, a configuration with Q 7^ gives 
positive contributions of the system size to the second 
and the third terms in the action (|10() . Thus we remove 
the first term from the action (|10|l . assuming that con- 
figurations with (5 7^ are negligible. 

The momentum cutoff A is determined by keeping the 
number of degrees of freedom for m. The variable m is 
originally defined for each hexagonal block with sides of 
length -\/3a (dash-dotted line in Fig. ^ and with area 
uo = (9-\/3/2)a^. Then, in the continuum limit, the mo- 
mentum cutoff A is given as aA = (2/3)(27r)^/^3~^/^ from 

{2TTf/L0 = 7rA2. 
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FIG. 2: Phase diagram in the space of distortion parameter 
5 and frustration parameter a (=J2/Ji). Phase boundaries 
1, and 



for 



|, and also the classical limit of the AF 
phase (5'=oo) are shown. For each 5", the region below the 
boundary is an AF order phase, while that above is a disor- 
dered spin-gap phase. Possible positions for InCu2/3Vi/303 
and Na3Cu2Sb06 are shown by diamond and triangle symbols 
respectively. 



We introduce rescaled dimensionless coordinates, a;o = 
Awo(c+c_)4r, xi — A(c_/c+)4x and X2 = A{c+/c^)^y 
with vo ~ SSaJi/ V2. The action 110|l with neglecting 
the Aq term is then rewritten in a standard NLSM form: 



1 

AeS = 



d^a 



dm 

dXa 



(13) 



with coupling constant go = 3^/2aAS^^{c+C- )~^ . 

Applying the one-loop RG analysis by Chakravarty et 
al. to Eq. H13|l . a quantum phase transition from the 
Neel phase with AF order to a disordered phase takes 
place at = Att. Rewriting this, the critical value ac for 
the phase boundary in the d-a space is given by 



1 



- ^6^(1-6)^ + X 



(14) 



with Xg — {^/3ttS^)^^ . Parameter As is a measure of 
quantum character for the spin system. In the classical 
limit of As = (5 = 00), Eq. (|14|l reduces to ac = 
i(l-(5) for (5 > and = ^{1 + 5-26^) for S < 0; the 
present theory is applicable below the classical line. The 
phase boundaries p4(l for 5 = i, 5* = 1 and S = ^ are 
shown in Fig. |21 For each S, the region above the solid 
line represented by the equation is a disordered phase 
with spin-gap and the region below is an AF ordered 
phase. Possible positions for materials InCu2/3Vi/303 
and Na3T2Sb06 are also shown. 
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In the case of no distortion (5 = 0), Eq. (|14() becomes 
simply a, = i(l - As); e. g. a, = 0.044 {S = i), 0.136 
{S = 1), and 0.153 {S = |). For S = ^, the disordered 
phase has a range of 74% against that of the classical 
AF order on the a-axis (Fig.[2J). The actual portion may 
be fairly reduced from this, since the NLSM method in- 
eludes, ;, g., a continuum approximation 

Whether 

or not a disordered phase survives for ^ = will be con- 
sidered in comparison with the square-lattice case. The 
portion of the disordered range in the J1-J2 model on a 
square lattice is (0.5-0.18)/0.5 x 100 = 64 % ||]. A disor- 
dered phase exists more plausibly for a honeycomb lattice 
than for a square lattice as was physically expected. For 
S' > 1, the disordered phase in Fig.|21is very narrow, and 
so its actual existence does not seem to be expected. 

In the limit of no frustration (a = 0), Eq. (|14|l is re- 
duced to 26^ - 3(5^ + 1 - A| 0. The critical value of 
distortion for S — ^ is 6 — 0.47. The disordered state 
on the (5-axis of Fig. |21 is explained by the tendency of 
forming singlet pairs at strong J{-bonds. 

Summarizing, we formulated an NLSM method for a 
spin system on a honeycomb lattice in the case that the 
exchange interaction includes both frustration by second- 
neighbor AF interactions and dimer-like distortion in 
first-neighbor AF interactions. Applying the RG analysis 
by Chakravarty et al. 0, we have a ground-state phase 
diagram in a parameter space. A disordered spin-gap 



phase is continuously extended from a strong frustration 
regime to a strong distortion regime, showing that the 
disordered ground states in both the regimes are essen- 
tially the same. In the case oi S — ^ , the spin-gap phase 
for the spin system on a honeycomb lattice is larger than 
that for the J1-J2 model on a square lattice, suggesting 
that a disordered state likely exists. 

Material NasTaSbOe with T=Cu {S = i) has a strong 
dimer-like distortion and a spin-gap unlike the other sim- 
ilar materials with T=Ni {S = 1) and T=Co (5" = |) 0. 
The distortion may be induced because it is in the spin- 
gap phase. In fact, within the spin-gap phase, the en- 
ergy of spins decreases and the lattice energy increases, 
as the dimer-like distortion increases; the balance deter- 
mines an actual distortion. On the other hand, within 
the AF phase, the dimer-like distortion does not reduce 
the spin energy, being consistent with little distortion in 
materials with T=Ni and T=Co. To find a material with 
small distortion in a spin-gap phase, it is reasonable to 
try S = ^ materials with strong elastic constants as well 
as appropriate second-neighbor AF interactions on hon- 
eycomb lattices. 
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